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ABSTRACT 

We have studied the dependence of the AGN nuclear radio (1.4 GHz) luminosity on 
both the AGN 2-10 keV X-ray and the host-galaxy K-band luminosity. A complete 
sample of 1268 X-ray selected AGN (both type 1 and type 2) has been used, which 
is the largest catalogue of AGN belonging to statistically well defined samples where 
radio, X and K band information exists. At variance with previous studies, radio upper 
limits have been statistically taken into account using a Bayesian Maximum Likelihood 
fitting method. It resulted that a good fit is obtained assuming a plane in the 3D Lr- 
Lx-Lk space, namely logL R = ^logL x + £ K log L K + £ , having a ~1 dex wide (la) 
spread in radio luminosity. As already shown by La Franca, Melini & Fiore (2010), no 



evidence of bimodality in the radio luminosity distribution was found and therefore 
any definition of radio loudness in AGN is arbitrary. Using scaling relations between 
the BH mass and the host galaxy K-band luminosity, we have also derived a new 
estimate of the BH fundamental plane (in the L^chz-^x-^bh space). Our analysis 
shows that previous measures of the BH fundamental plane are biased by ~0.8 dex 
in favor of the most luminous radio sources. Therefore, many AGN studies, where 
the BH fundamental plane is used to investigate how AGN regulate their radiative 
and mechanical luminosity as a function of the accretion rate, or many AGN/galaxy 
co-evolution models, where radio-feedback is computed using the AGN fundamental 
plane, should revise their conclusions. 

Key words: galaxies: active - radio continuum: galaxies - X-rays: galaxies - methods: 
statistical. 



1 INTRODUCTION 

In the last years, in many galaxy formation models, Ac- 
tive Galactic Nuclei (AGN) have been considered related 
to mechanisms capable to switch off the star formation in 
the most massive galaxies, thus reproducing both the ob- 
served shape of the galaxy luminosity function and the red, 
early type, passive evolving nature of local massive galax- 
ies. It is expected that the AGN and galaxy evolutions 
are closely connected to each other (the AGN/galaxy co- 
evolution) through feedback processes coupling both the star 
formation and the black hole (BH) accretion rate histories. 
Some of these models assume that the AGN feedback into 
the host galaxy is due to the kinetic energy released by the 
radio jets and is, therefore, dependent on the AGN radio lu- 
minosity ([Croton et aT||2006| [Cattaneo et al.||2006| parulli 



a kinetic luminosity function provides the adequate amount 



et a!|2008| ). It has been, indeed, already demonstrated that 
the conversion of the AGN radio luminosity function into 



of energy 


Best et al.|2006 Merloni & Heinz|2008 


Shankar 


et al. 2008 


Kording, Jester & Fender 2008 Cattaneo & Best 


2009 Smolcic et al.|2009 La Franca, Melini & Fiore|20 10). 



In this context, in order to build up more realistic AG- 
N/galaxy co-evolutionary models, it is very useful to mea- 
sure the dependence of the AGN core radio luminosity on 
other physical quantities, related to the AGN/galaxy evolu- 
tionary status, such as the BH and galaxy star masses and 
their time derivatives: accretion and star formation rates. 
These quantities can be, either directly or indirectly, mea- 
sured. 

A good estimate of galaxy star masses can be obtained 
by spectral energy distribution (SED) analyses in the optical 
and near infrared (NIR) domains (see e.g. Merloni et al. 
2010||Pozzi et al.|2012| >, or, still satisfactorily, from NIR (e.g. 
K-band) luminosity measures, as the mass-to-light ratio in 



the K band has a la scatter of 0.1 dex (Madau, Pozzetti & 
Dickinson|[l998} |Bell et~aT|2003 1. 
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BH masses can be estimated using reverberation map- 
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ping techniques or measuring the width of broad emission 
lines observed in optical and NIR spectra (the single epoch 
method; see e.g. Vestergaard 2002 1. Less direct estimates 



method as used by La Franca, Melini & Fiore (20101, we 



are obtained using the scaling relations observed between 
the BH mass and the bulge or spheroid mass, or between 
the BH mass and the bulge luminosity of the host galaxies 



(eg. |Dressler||l989 


Kormendy & McClure||1993 Kormendy 


& Richstone 1995 


Magorrian et al. 1998 1. In this framework, 



even the total host galaxy K-band luminosity, if converted 
into the bulge luminosity, can be used as a good proxy of 
the BH mass (see e.g. Fiore et al.|2012 l. 

The accretion rate, m, is related to the AGN hard (>2 
keV) X-ray luminosity, Lx, via the knowledge of the X-ray 
bolometric correction, Kx, and the efficiency, e, of conver- 
sion of mass accretion into radiation, 



Lx = 



(1) 



Lboi _ erne 
~Kx~ ~ (1 - e)K x ' 

where Lboi is the bolometric luminosity and typical values 



for e are about 0.1 ( jMarconi et al.|2004 Vasudevan & Fabian 
20091. 



The AGN radio luminosity distribution and its rela- 
tionship to either the optical or X-ray luminosity has been 
studied by many authors. Some of these studies discussed 
the AGN radio luminosity in terms of a bimodal distribu- 
tion where two separate populations of radio loud and radio 
quiet objects exist (e.g.|Kellerman n et al.|198 9; Mille rTF^ea-l 
cock & Mead 1990). Many other studies have alternatively 
measured a relationship between the radio and X-ray lumi- 
nosities (e.g. Brinkmann et ah] 2000| |Terashima fc Wilson 



2003} |Panessa et al.||2007| |Bianchi et al.||2009| |Singal et al 



2011 1. However almost all these studies are based on incom- 
plete samples due to the lack of deep radio observations. 
AGN samples selected in other bands (typically optical or 
X-ray) are therefore not fully detected in the radio band. 

However, in order to properly study the AGN radio lu- 
minosity properties and their relation to X-ray and optical 
luminosities, it is necessary to use fully (or almost) radio 
detected and complete AGN samples and when needed, to 
properly take into account the radio upper limits. More re- 
cently, using deep radio observations, it has been shown that 
both the AGN radio/optical and the radio/X-ray luminosity 



ratios span continuously more than 5 decades (Best et al. 
20051 l La Franca, Melini fc Fiore|[20Tol |Singal et al.||2011| 



Balokovic et al.|2012 i, without evidences of a bi- modal dis- 
tributions, and it is, therefore, inaccurate to deal with the 
AGN radio properties in terms of two separate populations 
of radio loud and radio quiet objects. La Franca, Melini & 
Fiore ( 2010 1 used a sample of about 1600 hard X-ray (mostly 
2-10 keV) selected AGN to measure (taking also into ac- 
count the presence of censored radio data) the probability 
distribution function (PDF) of the ratio between the nu- 
clear radio (1.4 GHz) and the X-ray luminosity Kx = log[v 
L„(1.4 GHz)/L x (2-10 keV)] (see |Terashima fc Wilson|2003[ 
for a discussion on the difference between the radio to op- 
tical and the radio to X-ray ratio distributions in AGN). 
The probability distribution function of Kx was functionally 
fitted as dependent on the X-ray luminosity and redshift, 
P(Rx\Lx , z). The measure of the probability distribution 
function of Kx eventually allowed to compute the AGN ki- 
netic luminosity function and the kinetic energy density. 
In this paper, using the same sample and a similar 



describe the measure of the dependence of the AGN radio 
core luminosity, Lr, distribution, on both the X-ray (2-10 
keV) luminosity, Lx , and the host galaxy (AGN subtracted) 
K-band luminosity, Lk- This measurement, in the context 
of the AGN/galaxy coevolution scenario (see above discus- 
sion), is vey useful in order to relate the kinetic (radio) feed- 
back to the accretion rate (Lx) and the galaxy assembled 
star mass (Lx). In order to accurately take into account 
the presence of censored data in the radio band, an ad hoc 
Bayesian Maximum Likelihood (ML) method and a three 
dimensional Kolmogorov Smirnov test have been developed. 

Many authors have observed the existence of an anal- 
ogous relationship between the radio luminosity, the X-ray 
luminosity, and black hole mass (M), the BH fundamental 
plane, namely logL/i = £rx logLx + ^RM^ogM + constant 



(see Merloni, Heinz & di Mattco 2003 Falcke, Kording fc 
Markoff|2004||Gultekin et al.|2009a| >. The measure of such a 
relation is very useful in order to discriminate among several 
theoretical models of jet production in AGN, as it suggests 
that BH regulate their radiative and mechanical luminosity 
in the same way at any given accretion rate scaled to Ed- 
dington (|Falcke fc Biermann||1995| |Heinz fc Sunyaev|[2003 



Churazov et al.|2005[ |Laor fc Behar||2008 1 

As relations have been observed between the BH mass 
and the K-band bulge luminosity (see discussion above), 
in section 6 we convert our measure of the relation in the 
logL_R-logLx-logL_ft: space, into a relation into the logL_R- 
logLx-logAf space, and discuss how much important is to 
properly take into account the presence of radio upper limits 
in order to measure the BH fundamental plane. 

Unless otherwise stated, all quoted errors are at the 
68% confidence level. We assume Ho=70 Km s _1 Mpc -1 , 
Q m =0.3 and Ov=0.7. 



2 THE DATA 



In our analysis we used the same data-set used by |La Franca,| 
Melini & Fiore (20101, where radio (1.4 GHz) observations 
(either detections or upper limits) were collected for 1641 
AGN (both type 1, AGN1, and type 2, AGN2, i.e. showing or 
not the broad line region in their optical spectra) belonging 
to complete (i.e. with almost all redshift and N^j measures 
available) hard X-ray (>2 keV; mostly 2-10 keV) selected 
AGN samples, with unabsorbed 2-10 keV luminosities larger 
than 10 42 erg/fQ 

As the goal was to use the radio luminosity in order to 
estimate the kinetic luminosity, |La Franca, Melini fc Fiore| 
( |2010[ ) measured a radio emission which was as much as pos- 
sible causally linked (contemporary) to the observed X-ray 
activity (accretion). Radio fluxes were measured in a region 
as close as possible to the AGN, therefore minimizing the 
contribution of objects like the radio lobes in FRII sources 



(Fanaroff & Riley 19741. For this reason they built up a large 
data-set of X-ray selected AGN (where redshift and N^j col- 
umn densities estimates were available) observed at 1.4 GHz 



1 Throughout this work we assumed that all the X-ray sources 
having 2-10 keV unabsorbed luminosities larger than 10 42 erg/s 
are AGN. See e.g. |Ranalli, Comastri fc Setti| pOOty for a study 
of the typical X-ray luminosities of star forming galaxies. 
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with a ~1" typical spatial resolution (1" corresponds, at 
maximum, to about 8 Kpc at z~2). The cross correlation 
of the X-ray and radio catalogues was carried out inside a 
region with 5" of radius (almost less than, or equal to, the 
size of the central part of a galaxy like ours), following a 
maximum likelihood algorithm as described by Sutherland 
fc Saunders] < fl992| > and |Ciliegi eTaL] p003| . The off-sets be- 
tween the X-ray and radio positions of the whole sample 
resulted to have a root mean square (rms) of 1.4" (similar 
to the typical values obtained in X-ray to optical cross- 
correlations; e.g. Cocchia et al.(2007 l. 

In order to measure the K-band galaxy luminosity, the 



1641 AGN from La Franca, Melini & Fiore (20101 have been 



cross-correlated with already existing K-band photometric 
catalogues as explained below. 
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4 


COSMOS 
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575 


121 


CLANS 


139 


125 


91 


52 


ELAIS 


421 


363 


283 


37 


CDF-S 


94 


86 


85 


12 


CDF-N 


162 


161 


158 


40 


Total 


1641 


1517 


1268 


293 



Table 1. Samples breakdown 



2.1 The local Sample: SWIFT and Grossan 

At the lowest redshift we have used a sample of 33 AGN 



belonging to the 22 month SWIFT catalogue ( Tueller et al 



20081. These AGNs have been detected at high galactic lat- 
itude (|fe| > 15deg) in the 14-195 keV band with fluxes 
brighter than 10 -11 erg s _1 cm -2 ; all objects have Nh col- 
umn density and optical spectroscopic classification avail- 
able. The radio luminosity at 1.4 GHz has been derived us- 
ing the Faint Image of the Radio Sky at Twenty cm (FIRST) 



Very Large Array (VLA) survey ( jBecker, White fc Helfarid| 
1995 1. In the case of no radio detection, a 5a upper limit of 



0.75 mjy was adopted. The correlation between catalogues 
was made through the likelihood ratio technique ( |Suther-| 
land fe Saunders||1992| |Ciliegi et~aL|2003[ ). K s magnitudes 
have been associated to the sources using the 2MASS (Two 
Micron All Sky Survey) catalogue which has a ~14.3 mag 
completeness limit magnitude (Skrutskie et al. 2006). The 



cross-correlation with the K-band data has been carried out 
by comparing the positions of the the 2MASS sources on the 
K-band image of the AGN counterpart. 

To enlarge the local sample, |La Franca, Melini fc Fiore] 
( |2010p used the hard X-ray selected AGN catalogue de- 
tected by the HEAO-1 mission (2-10 keV fluxes b r ighter 
than 2 x 10 -11 erg s _1 cm -2 ) described by |Grossan| ( |1992[ ) 
and revised by |Brusadin| ( [20031 ). As for the SWIFT sample 
these sources have been cross correlated with the FIRST 
and 2MASS radio and K-band catalogues, respectively. 

In summary, the local sample contains 43 X-ray sources 
all having a Ks band detection. 



2.2 HBSS 

We have used the 32 AGN selected by the Hard Bright Sen- 
sitivity Survey of the XMM-Newton satellite, carried out 
at 4.5-7.5 keV fluxes brighter than 7xl0 -14 erg s _1 cm -2 
(Delia Ceca et al. 2008). Ks detections and upper limits 



were obtained for 19 and 13 sources, respectively, from the 
2MASS catalogue. The cross-correlation was carried using 
the same technique used for the local sample. Those sources 
missing a Kg detection were eventually excluded from the 
analysis. 



2.3 The ASCA surveys: AMSS and ALSS 

Two samples come from observations of the ASCA satel- 
lite: the ASCA Medium Sensitivity Survey (AMSS; |Akiyama| 
et~aLl[2003t , which is com posed of 43 AGN, and the ASCA 
Large Sky Survey (LASS; Ueda et al.|fl999 ) , which is com- 



posed of 30 AGN. For both of these catalogues we used the 



Ks photometric measures by Watanabe et al. (20041. Only 



one source (belonging to the AMSS) misses the Ks detec- 
tion. 



2.4 COSMOS 

COSMOS is the biggest catalogue used in this work. The cat- 
alogue by Cappelluti et al.| (|2009|) of XMM-Newton sources 

-2 



with 2-10 keV fluxes brighter than ~3xl0 15 erg s 1 cm 
was used. K-band magnitudes for 648 out of 677 sources were 



taken from Brusa et al. (2010 and private communication) 



2.5 CLANS 

K-band photometry for 125 out of 139 AGN belonging to 
the Chandra Lockman Area North Survey were obtained 



Morrison (20081 



from Trouille et al. (20081. The radio data comes Owen & 



2.6 ELAIS-S1 

In the field SI of the European Large Area ISO Survey 
(ELAIS-S1) we used the catalogue of X-ray sources detected 
by XMM-Newton by |Puccetti et al.| ( [2006| which reaches a 2- 
10 keV flux limit of 2x 10 -1 ° erg s _1 cm -2 . Radio data were 
taken from Middclberg et al. ( 2008 1, while spectroscopic and 



photometric identifications were taken from|La Franca et al 



20041 ;|Berta et al. (20061; Feruglio et al.|p008l and Sacchi 



et al.| ( |2009[ ). In |Feruglio et al.| ( |2008[ ) K s detections for 363 
objects out of the 421 sources used by |La Franca, Melini fc| 
Fiore (20101 are available. 



2.7 Deep samples: CDF South and North 

The deepest X-ray catalogues used in this work are those 
available in the Chandra Deep Field South (CDFS) and 
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North (CDFN). These are sub-samples of the GOOD-S and 
GOOD-N multi-wavelength surveys with 2-10 keV flux lim- 
its of 2.6 xlO" 16 ergs" 1 cm" 2 and 1.4xl0" le ergs" 1 cm" 2 , 
respectively. 

In the CDFS we used 94 sources from the catalogue of 



Alexander et al. (20031 and identified by Brusa et al. (2010) 



These sources were correlated with radio data taken from 



Miller et al. (2008). Ks band photometry was obtained us- 



ing the deepest multi-wavelength catalogue (FIREWORKS) 
provided by Wuyts et al. ( 2008 1 that reaches Ks —22.5 mag. 
The Ks band catalogue was cross correlated with the opti- 
cal catalogue using the likelihood ratio technique. Ks band 
counterparts for 86 out of 94 sources were found. 

In the CDFN we used the X-ray catalogue from Alexan- 



der et al. (2003) with the identifications from Trouille et al. 
(2008). The sample consists of 162 extragalactic sources for 
which radio informations were obtained from Big gs fc Ivison| 
(20061. iTrouille et al.l (120081) provided K s band detections 



for all but one of the sources. 

In Table [l] the breakdown of all samples used is re- 
ported. In column 1) it is shown the original number, N, 
of the sources contained in the samples used by |La Franca, | 



Melini & Fiore (20101, while in column 2) it is shown the 



number of sources, Njr, having a K-band detection. 



3 K-BAND AGN AND GALAXY 
LUMINOSITIES 

The K-band absolute magnitudes, Mk, have been computed 
by applying an empirical K-correction. We used the formula 



m x +5- 51ogd,(z) + 2.5(1 + a)log(l + z), (2) 



where di(z) is the luminosity distance and the K-correction 
is represented by the a parameter. We used a — —0.86 af- 
ter a comparison of our data with the COSMOS catalogue 
where absolute K-band magnitudes have been accurately 



computed through SED fitting techniques (Bongiorno et al 



2012 and references therein). 



In order to estimate the host galaxy K-band luminosi- 
ties (i.e. the stellar component) we subtracted the AGN con- 
tribution from the measured total luminosities. For this pur- 
pose we used the nuclear (AGN only) infrared SEDs, nor- 
malized to the hard X-ray (2-10 keV) intrinsic luminosity 
and averaged within bins of absorbing Njj as published by 



Silva, Maiolino & Granato (20041. These AGN SEDs were 



obtained trough the interpolation, using updated models 



from Granato & Danese (19941, of the nuclear IR data of 



AGN taken from the Maiolino & Rieke ( 1995 1 sample. The 



accuracy of our method has been tested by comparing our 



estimates with those obtained by Merloni et al. ( 2010 1 and 
|Bongiorno et al. (20121 on the COSMOS catalogue using 
SED decomposition fitting techniques. As shown in Figure 
[I] our estimates, although less accurate, are in good agree- 



ment with those obtained by Merloni et al. (2010 \ and Bon 
|giorno et al. (2012). The average difference results to be 
logL /f (COSMOS)-logL K (us)=-0.07 (0.05) dex, with a la 
spread of 0.30 (0.18) dex for AGN1 (AGN2). 

In some cases we obtained that the expected AGN K- 
band luminosity was very close to (or even larger than) the 
total measured (AGN + host galaxy) luminosity. With the 
SED decomposition fitting techniques used by Merloni et al. 



100 



mean = 0.054 ct = 0.18 
mean=-0.07 <j = 0.30 




-0.5 0.5 

log(L k (Bongiorno)/L k (us)) 



Figure 1. Histogram of the logarithmic differences between our 
estimates of the K-band luminosities of the host galaxies (after the 
AGN component subtraction; see text) in the COSMOS sample 
and the SED fitting measures from Bongiorno et al. (2012). AGN1 
and AGN2 are shown by blue and red lines, respectively. 



(2010) and Bongiorno et al. (20121 no object could result to 



have an AGN luminosity larger than the total one (see also 
|Pozzi et al.][2007l [20T2] ). |Bongiorno et al.| poHf conserva- 
tively decided that if the galaxy component were smaller 
than 10% of the total one, only an upper limit, correspond- 
ing to 10% of the total luminosity, could be assigned. Fol- 
lowing this approach, we decided to adopt a more conserva- 
tive assumption, and excluded from our analysis those 249 
sources where the galaxy component resulted to be smaller 
than 20% of the total on<£] 

In Table[l]we report in column 3) the number, Nk-cl, 
of sources where it was possible to estimate the galaxy K- 
band luminosity. 



4 THE WHOLE SAMPLE 

In summary our data-set is composed by 9 X-ray selected 
AGN samples which contain a total of 1268 sources for which 
we were able to estimate the K-band stellar component lu- 
minosities of the host galaxy, Lk- For all these AGN, column 
densities and de-absorbed 2-10 keV luminosities, Lx, mea- 
sures are available. The radio data allowed to measure on 
293 sources the "nuclear" 1.4 GHz luminosity, L_h, while for 
the remaining sources 5a radio upper limits (Table [l] col- 
umn 4) are available. In total, this is the largest catalogue of 
AGN (both of type 1 and 2) belonging to statistically well 



2 As described in section 5, our 3D ML fitting method is able to 
deal with upper limits on one physical quantity only (the radio 
luminosity in our case). Therefore all sources where an upper limit 
on their K-band luminosity was available, were excluded from our 
analysis (this happened to all the 10 sources of the GROSSAN 
sample). 
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Figure 2. 2-10 keV de-absorbed luminosity, Lx, of the total sam- 
ple as a function of redshift . 
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Figure 3. Distribution of the sample onto the 2D Ljf-Lx, ^k' 
hji and Lx-Lfl planes, where L^* is the total (galaxy + AGN) 
K-band luminosity. Squares represent the radio detected sources 
while the sources with a radio upper limit are shown by either 
grey arrows or grey squares (upper left panel). 



defined samples where radio, X and K band information ex- 
ist. The distribution, in the Lx-z plane of the whole AGN 
sample is shown in Figure [2] while in Figure [3] and Figure 
[4] we show the 3D Lk-Lx-L_r distribution, projected on to 
the three 2D Lk-Ljc, Lk-L_r and Lx-L/j planes, where Lk 
is shown before (L^*), and after (Lk), the subtraction of 
the AGN component in the K-band, respectively. 



Figure 4. As in Figure [3] with Ljf 
tracted K-band galaxy luminosity. 



representing the AGN sub- 



5 THE PLANE FITTING 

5.1 Maximum Likelihood fitting method 

We have used a maximum likelihood fitting technique, with 
a Bayesian approach, in order to derive the probability dis- 
tribution function of the AGN radio luminosity, L_r, as a 
function of Lx and the K band stellar component lumi- 
nosity, Lk, P(Lr\Lx, Lk)- The maximum likelihood fitting 
method doesn't need to use binning (as, e.g., it happens 
when using the \ 2 fitting method), and therefore the results 
do not depend on the arrangement of the binning pattern. In 
the usual maximum likelihood fitting of luminosity function 
distributions p(z,L) of extragalactic sources 



(eg. 



Marshall 



et al.|[T983| l the best fit solution is obtained by minimizing 
the quantity S = —21nC (where C is the likelihood, and the 
S function follows the x 2 statistic and therefore allows to 
estimates the confidence interval of the best fit parameters; 
Lampton, Margon & Bowycr 1976). The natural logarithm 
of the likelihood function is computed as follows: 



InC = ^2 ln(p(zi, Li)) — J 



p{z,L)D.(z,L)'^dLdz, (3) 



where the sum is made over all the i observed sources 
and Q,(z,L) is the sky coverage as function of the luminos- 
ity L and redshift z. The first term is proportional to the 
combined probability (of independent events) of observing 
all the i sources, each having redshift Zi and luminosity Li, 
while the second term corresponds to the total number of 
expected sources in the sample and is, therefore, also used 
to constrain the normalization of the luminosity function. 

In our case we have devised a new S function to be 
minimized which, following the same statistic principles of 
the maximum likelihood method, is able to measure the 
conditional probability distribution function P(Lr\Lx, Lk) 
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(where j P(Lr\Lx , Lk)(ILr = 1) of observing an Lr radio 
luminosity in an object having Lx and Lk luminosities. Our 
method has the advantage (in comparison to other three- 
dimensional fitting techniques) to be also able to take into 
account the occurrence of upper limits in one of the three 
dimensions. Indeed, as already discussed, in all samples the 
radio observations are not deep enough do detect all the 
sources (see Table [T| and therefore upper limits need to be 
considered in order to derive the true AGN radio luminos- 
ity distribution (see |Plotkin et al. 2012 for an analogous 
Bayesian approach on this topic). For these reasons the nat- 
ural logarithm of the likelihood function has been computed 
as follows: 

lnC = ^lnP{L Ri \Lx i ,L Ki ) - 



J P{L R \L X] ,L K] )dL R . 

j L R >L lim 



(4) 



The first sum is computed for all the i radio-detected 
sources and, as in eq. [3j is proportional to the combined 
probability of observing the radio luminosities L Ri of all 
the i detected sources, while the second term is the sum of 
the probability of radio detecting each j observed (either 
radio detected or not) AGN, with a radio luminosity larger 
than its radio detection limit Ln mj . In this case, analogously 
to the classical maximum likelihood method (eq. [jjjl, this 
second term corresponds to the expected total number of 
radio detected sources. 



5.2 The 3D Kolmogorov-Smirnov test 

Although the maximum likelihood technique is very pow- 
erful in finding the parameters of the best fit solution and 
their uncertainties, it doesn't allow to quantify how good 
the solution is. We have then devised a three-dimensional 
Kolmogorv-Smirnov (3D-KS) test able to measure the prob- 
ability of observing the 3D distribution (in the Lr—Lx—Lk 
space) of our sample of radio detected sources, under the null 
hypothesis that the data are drawn from the best fit model 
distribution. 

The K-S test is a standard statistical test for deciding 
whether a set of data is consistent with a given probabil- 
ity distribution. In one dimension the K-S statistic is the 
maximum difference D between the cumulative distribution 
functions of the data and the model (or another sample). 
What makes the K-S statistic useful is that its distribution 
(in the case of the null hypothesis that the data are drawn 
from the same distribution) can be calculated giving the 
significance of any nonzero value of D. 

In order to make a 3D-KS test we have followed the 
examples of the generalization of the K-S test to two- 



dimensional distributions by Peacock (1983) and Fasano & 



Franceschini ( 1987 1. The maximum difference D statistic has 



been computed by measuring the difference among the frac- 
tion of observed and expected (by the model) sources in 
each of the eight octants defined at the 3D positions of each 
radio detected source. In order to compute the number of 
expected sources, 2000 Monte Carlo simulations have been 



used. Although the 3D-KS test is carried out on the radio 
detected sample only, it properly takes into accounts the ef- 
fects of the radio upper limits. Indeed, the simulations have 
been carried by extracting a radio luminosity (or not) for 
each observed source by taking into account its radio detec- 
tion limits an the model conditional probability distribution 
function P(Lr\Lx , Lk)- The same 2000 simulations have 
been used to compute the probability (significance) to ob- 
serve the measured maximum D statistic in the case of the 
null hypothesis that the data are drawn from the same dis- 
tribution. 



5.3 The fits 

We have tried to fit the data assuming that the AGN radio 
luminosity, Lr, depends on average (i.e. with a spread in the 
radio luminosity axis) linearly from both the logarithmic X- 
ray and K-band (host galaxy) luminosities, drawing a plane 
in the 3D \ogL R-\ogL x-\ogL k space. As already discussed 
in the introduction, this assumption is suggested by the ob- 
servations, done by many authors, of the existence of a BH 
fundamental plane: a similar relationship between the radio 
luminosity, the X-ray luminosity, and black hole mass (A/), 
namely logLfl = ^Rx^ogLx +^RM^ogM + constant (see Mer- 



loni, Heinz fc di Matteo|[2003} [Gultekin et al. 2009a). The 
AGN radio luminosity dependence on logLx and logL^ was 
modeled by the following relationship: 



logL H = £xlogLx,44 + £A"logL K ,43 + £o, 



where L r is the luminosity of the peak (mode) of the 
probability distribution function of the spread in erg s _1 
units, and the luminosities have been normalized to Lx = 
10 44 erg s _1 Lx,44 and L K = 10 43 erg s _1 Lk,43- 

This spread was first modeled by a Gaussian function, 
G(X\0, a), where X — Io^Lr/Lr), centered on X = with 
standard deviation a. The best fit parameters are shown in 
Table [5] (model 1). However, the 3D-KS statistic tells that 
there is only a 2% probability that the observed distribution 
is drawn from the model. The probability distribution func- 
tion of X is shown in Fig. [5] The data have been plotted 
comparing in each bin of X the number of observed (N bs) 
and expected (N exp ; by the model) sources. This method 
(N i,s vs. N exp method) reproduces the observations and 
consequently properly takes into account both the radio de- 
tections and the upper limits (see e.g. La Franca et al.|1994 



|La Franca fc CristiariTl[T997l |Matute et al.||2006| for simi- 
lar applications). It is worth to note that this fit, although 
not satisfactory, gives and indication that the PDF of Lr 
is quite large: it results a ~ 1 dex. i.e. 68% of the cases 
are included into a 2 dex wide distribution. Similar results 



have been found by Merloni, Heinz & di Matteo ( 2003 1 and 
|Giiltekin et~aT] ( |2009a[ H 



As shown in figur e |5| (a nd already observed by La 



Franca, Melini & Fiore 20101, the data show an excess of 



high radio luminosity sources if compared to a symmetri- 
cal distribution. In order to better take into account this 
excess we have assumed an a-symmetrical double Gaussian 
distribution with two a values: 01 and a u for radio luminosi- 
ties values below (lower) and above (upper) the value of the 
radio luminosity, L r , of the peak of the probability distri- 
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Model 


Function 


5x 




So 


Off-set 






a 


k 


S 




1 

2 
3 


Gaussian 
Double Gaussian 
Lorentz 


0.313 
0.373 
0.379 


0.683 
0.582 
0.705 


38.684 
38.513 
39.374 


0.00 
0.41 
-0.43 


0.974 
0.497 
1.141 


1.02 
0.30 






1583 
1571 
1515 


2% 
3% 
26% 


4 
1 


Gaussian + exp 
68% conf. errors 


0.387 

+0.031 
-0.059 


0.632 

+0.066 
-0.053 


38.937 

+0.050 
-0.032 


0.09 


0.582 

+0.028 
-0.038 




0.117 

+0.087 
-0.021 


1.66 

+0.24 
-0.07 


1504 


35% 



1) Corresponding to <r in those models with a single spread parameter. 



Table 2. Best fit solutions 



-3 




Gaussian+exp 




-1 1 

X=log(L R /L R ) 



Figure 5. Probability distribution functions of X = 
log(L r I L r) , where L r is the luminosity of the peak {mode), 
whose position is shown by a vertical continuous orange line. The 
four distributions correspond to the four best fits reported in Ta- 
ble^ The blue vertical lines show the values of the means of the 
distributions, whose off-set from the peak luminosity is reported 
m Table [5] and which are used in the next figures to represent 
the best fit solutions. 



log L R [erg s 1 
44i 




log L x [ergs" 1 ] 42 47 



log L K [erg s x ] 



Figure 6. 3D distribution in the L^j-Lx-Lif space. Filled circles 
represent radio detections and faint open triangles represent radio 
upper limits. The plane of the distribution of the mean Ljj as a 
function of Lx and L/f , according to the best fit solution 4 (Table 
|2}, is shown. 



the fit with the data) should, indeed, be carried out using 
the position of the mean of the distributions. The offsets 
between the mean and the mode (off-set = mean - mode) of 
the PDF are also listed in column 4 of Table [2] 

Better results are obtained if the spread is modeled, as 
proposed by La Franca, Melini & Fiore (20101, assuming a 



double Lorentzian function described by the parameters 07 
and a u : 



butiorj^] (i.e. for X values above or below zero) respectively. 
Indeed, the fit gives a larger spread, cr„ = 1.0, at X>0, than 
measured at X<0, where cr;= 0.5. However, even with this 
model, the 3D-KS test gives a not good enough, 3%, prob- 
ability. As the PDF is a-symmetrical, we show in Figure [5] 
both the locus of the peak and of the mean of the PDF. A 
better representation of the model (e.g. in figures comparing 



defined by eq. [5] 



P(X) = 4 



(X > 0), 



(6) 



-car 



where, in order to obtain a continuous function at X=0, 
A = \/ G u j oil an d the parameter TV is constrained by the 
probability normalization requirement: j P(X)dX = 1. In 
this case (see Figure|5]and model 3 in Table[2) a 26% 3D-KS 
probability is obtained. 
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44 



42 



JL 40 



M 
O 






SWIFT 

CLANS 



-2 2 

0.387 log L X44 + 0.632 log L K43 

Figure 7. Edge-on view of the 3D plane. The more complete 
samples, SWIFT and CLANS, are represented by red and green 
squares, respectively. The continuous line shows the locus of the 
mean of the radio luminosity probability distribution function of 
our best fit solution 4 (Table [2| . Radio upper limits are repre- 
sented by arrows. 



An even better solution is obtained if, for X > 0, we 
add to a Gaussian distribution, G(X |0, a) (as in model 1), 
an exponential function able to reproduce the high radio 
luminosity tail, 



P(X) 



aX 2 



bG(X\0,a) 
bG(X\0,a) 



(X<0) 
(X > 0), 



(7) 



where the a and b parameters are not independent, as they 
are constrained by the probability normalization require- 
ment 



aX 2 e~ kx dX. 



(8) 



The best fit solution (model 4 in Tableland see Fig- 
ure [5]l gives b=0.9499, which implies that the added expo- 
nential tail at high radio luminosities represents about 5% 
(1-0.9499) of the population (10% for XJsO where it is de- 
fined). In this case a 35% 3D-KS probability is obtained. 



As already shown by La Franca, Melini & Fiore (20101 the 



radio luminosity distribution of the AGN does not show any 
evidence of a bimodal distribution and therefore any defi- 
nition of radio loudness is arbitrary. Nonetheless, it should 
be observed that (as demonstrated by these fits) the radio 
luminosity PDF is not symmetrical, but skewed with a long 
tail at high radio luminosities. The introduction of this tail 
at X>0 allows to obtain a narrower (a ~ 0.6 dex) comple- 
mentary symmetrical Gaussian distribution. The position 
of the peak of the radio distribution of our best fit solution 
(number 4) is represented by the following equation: 



logLfl = 0.39t^logL x ,44 + 0.63±-{£logL*r,43 + 39.94±"g|, 

(9) 

Confidence regions of each parameter were obtained by 
minimizing the S function at a number of values around 



41 



40.5 



40 



39 



M 
o 



38.5 



38 



37.5 



"i r 



"i i 

SWIFT 

AMSS 

ALSS 

HBSS 

COSMOS 

CLANS 

ELAIS 

CDF-S 

CDF-N 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Completeness Fraction 



Figure 8. Mean radio luminosity of the radio detected sources 
of each sample as a function of the radio detection completeness. 



the best-fit solution, while leaving the other parameters free 
Lampton, Margon fc Bowyer||1976"l ). The 68% 



to float (see 

confidence regions quoted correspond to AS (=A^ 2 ) = 1. 
The best fit solution is also shown in Figure [6] where the 
3D distribution of the data is shown, and Figure [jj where 
the 2D edge on view of the plane is shown. This last fig- 
ure helps understanding how much important is to take into 
account the effects of using censored data. The fitting solu- 
tion seems, indeed, to be not a good representation of the 
distribution of the radio detections. This is because most of 
the X-ray samples have no radio observations able to de- 
tect all sources and then the radio detections are biased in 
favor of the most luminous radio sources. Indeed the mean 
radio luminosity of each of our samples decreases as a func- 
tion of the radio identification completeness fraction: the 
average radio luminosity of the samples with about 10% ra- 
dio identifications is of about logLj;=40.0 erg/s, while for 
the 90% complete samples the average luminosity is about 
2 dex lower (logL_H=38.2 erg s _1 ; see Figure [s|. This bias 
is properly taken into account with our ML fitting method 
which takes into account upper limits. The fitting solution, 
indeed, fairly well reproduces the distribution of the most 
radio complete samples, such as SWIFT and CLANS (see 
Figures [jj and ^ . 



6 THE AGN FUNDAMENTAL PLANE 

Many authors have observed the existence of a BH funda- 
mental plane relationship between the radio luminosity, the 
X-ray luminosity, and the BH mass (see the discussion in the 
Introduction). As a relation has been observed between the 
BH mass and the bulge luminosity, it is interesting to study 
whether the BH fundamental plane measures are compatible 
or not with our measured relationship in the logL_R-logLj>f- 
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estimates with those reported by Merloni et al. ( 2010 1, which 



log L R [erg s 1 
44i 




43 

log L x [ergs -1 



log L K [erg s *] 



Figure 9. Same as Figure [fj] but using only the, more radio com- 
plete, SWIFT (red circles) and CLANS (green circles) samples. 
Radio upper limits are represented by open triangles. The plane 
of the distribution of the mean Lr as a function of Lx and Lk, 
according to the best fit solution 4, is shown. 



logL/c plane (we remember that Lk is the galaxy, stellar 
component, luminosity). We have therefore estimated the 
BH masses using the calibrated black hole versus K-band 



bulge luminosity relation from Graham (20071 
log(M bh /M e ) = -0.37(M K ,buise + 24) + I 



.29. 



(10) 



The bulge to total luminosity ratio (B/T), or the bulge to 
disc ratio (B/D; where B/T = [1 + D/B]' 1 ) are function 
of the galaxy type (see |Dong fc De Rob ertis 2006} |Graham| 
& Worlcy 2008). Therefore, in order to properly derive the 
bulge luminosities from our measures of the total galaxy lu- 
minosities it would be necessary to know the morphological 
type of our galaxies. Unfortunately this information is barely 
available only for few tenths of galaxies belonging to the lo- 
cal sample (SWIFT), while for the higher redshift galaxies 
no information is available. Moreover, such a kind of rela- 
tionships have been calibrated only on local samples, while 
it is well known that the average galaxy mass and morphol- 
ogy changes with redshifts (e.g. smaller and more clumpy 
galaxies with increasing redshift; Mosleh et al. 20121, and 
our sample reaches z ~ 5. However, our aim is not to mea- 
sure the BH fundamental plane, but only to verify how much 
compatible our measure is with previous measures. Accord- 
ing to Graham fc Worley] ( |2008 1 , we have therefore assumed 
an average value of B/T = 1/4 in the K bandpass (see e.g. 
Fiore et al. 2012| for similar assumptions). Under this as- 
sumption, eq. |10| corresponds to the relation 



log(M 6fe /M ) = 0.9251ogLjr - 31.781, 



(11) 



where Lk is the total (bulge plus disk) galaxy luminosity, 
expressed in erg s _1 units. We have compared our BH mass 



have been obtained via virial based analysis of optical spec- 
tra of AGN1. On average, our estimates are larger by 0.14 
dex in solar masses with a spread of 0.43 dex, in solar mass 
units. This spread is compatible with the typical spreads 
(~0.3-0.4 dex) of the BH mass estimates based on both 
the virial and the scaling relations methods ( Giiltekin et al 



2009b I, whose uncertainties, in our comparison, should be 



both taken into account in the propagation of errors. 

According to eq. 1 1 1 1 our best fit solution 4 is trans- 
formed into the L_r-Lx-M_bh space in to the relation 

logL fl = 0.391ogL x + 0.681ogA/ BH + 16.61 (12) 



while Merloni, Heinz & di Matteo (2003) measure 



logL R = 0.601ogLx + 0.781ogA/ S H + 7.33, 



(13) 



where Lr here is the 5 GHz nuclear luminosity in units erg 
s _1 , with 1.4 GHz radio luminosities converted into 5 GHz 
vL v luminosities assuming a radio spectral index a — 0.7 
(where L„ oc v' -°; |Condon, Cotton fc Broderick||2002| |, L x 
the 2-10 keV nuclear X-ray luminosity in units of erg s _1 , 
and Mbh the black hole's mass in units M© (Merloni, Heinz 
fc di Matteo|2003 |. 



7 DISCUSSION AND CONCLUSIONS 



As expected (see Figure 10 1, our estimate of the BH funda- 
mental plane predicts lower radio luminosities if compared 
with t he previous measure by |Merloni, Heinz fc di Matteo| 



1.2003). The typical difference, computed at 10 Mq BH mass 



and 10 44 ergs^ 1 X -ray luminosity, is of ~0.8 dex. As already 
observed in the Lr-Lx-Lk space, this difference is due to 
the inclusion in our analysis of the contribution of the radio 
upper limits, and indeed our best fit solution reproduces well 
the distribution of the most radio complete samples such as 
SWIFT and CLANS (Figures 



11 



and 



12) 



It should be noted that the fundamental plane of |Mer-| 
loni, Heinz fc di Matteoj (|2003[) was constructed by including 



a sample of X-ray BH binaries. While it is interesting to see 
that the mass scaling is still broadly consistent with that 



of Merloni, Heinz & di Matteo ( 2003 1 even in our study 



of an AGN-only sample, the radio-X-ray coefficient is very 
different: ~0.4 in this study, while low-accretion rate (and 
low-luminosity) X-ray BH binaries show a tight radio-X- 
ray correlation with slope ~0.6. Interestingly enough, re- 
cent re-analysis of radio-X-rays correlations in X-ray bi- 
naries suggests the presence of a second, less radio lumi- 



nous branch (Gallo, Fender & Pooley 2003 Coriat et al 
2011} |Gallo, Miller fc Fender||2012[ ). In this framework, our 



study, could suggest that, also in AGN, a second less radio 
luminous population should be taken into account, which 
could corresponds to those objects with low values of the 
X — \og{Ln/ Lr) parameter (see Figure[5}. 

As already discussed in the Introduction, the measure of 
the dependence in AGN of Lr from other physical quantities 
such as Lx and the host galaxy K-band luminosity, Lk, is 
very useful in order to better include AGN in the galaxy 
evolution models where the AGN/galaxy feedback play a 
relevant role. 
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Figure 10. 3D distribution in the in Ljj-Lj-Mg^ space. The 
projection from the L^j-Lx-Lk space of the plane of the best fit 
solution 4 is shown in cyan. The fundamental plane from | Merloni, | 
|Heinz fc di Matteo| ( |2003| > is shown in gold. Radio upper limits 
are represented by open triangles. 
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Figure 11. Same as Figure [To] but using only the, more radio 
complete, SWIFT (red circles) and CLANS (green circles) sam- 
ples. Radio upper limits are represented by open triangles. 
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Figure 12. Edge on view of the BH fundamental plane as mea- 
sured by |Merloni, Heinz fc di Matteoj ( |2003} (gold dashed line) . 
The continuous cyan line shows our solution for a sample of AGN 
having average logM=8.5. The more complete samples, SWIFT 
and CLANS, are represented by red and green squares, respec- 
tively. Radio upper limits are represented by arrows. 



Our analysis has allowed to find a good analitical solu- 
tion represented by a plane in the 3D logL/j-logLx-logLx 
space once a wide (1<j~1 dex), a-symmetrical, spread in the 
radio luminosity axis is included. This result confirms the 
study of La Franca, Melini & Fiore (20101, who, studying 



the dependence of PDF of Lu from Lx and z, where able to 
model the 1 dex wide (la) spread in the AGN radio lumi- 
nosity distribution. This results show that a proper study of 
the correlation between different bands luminosities in AGN 
(or other sources) cannot be performed without taking into 
account censored data. 

A clear example, in this framework, is the measure of 
the BH fundamental plane in the 3D logLn-logLx-logM 
space. After converting, using scaling relations, our mea- 
sures of the host galaxy K-band luminosity into BH masses, 
our best fit solution corresponds to a BH fundamental plane 
which on average predicts 0.8 dex lower values for the AGN 
radio luminosities. 

It should be pointed out that, at variance with many 
similar statistical studies, our analysis is based on a compila- 
tion of complete, hard X-ray selected, AGN samples, where 
both AGN1 and AGN2 are included. Therefore, our results 
should better represent the behaviour of the whole AGN 
population. However it will be interesting to check the 3D 
correlation between Lr, Lx and Lk in complete, radio se- 
lected samples. As the Merloni, Heinz & di Matteo (20031 



sample was a hybrid sample without a clear selection crite- 
rion, it is plausible that parts of the differences found in the 
present work could be due to the specific selection criterion. 
In very general terms, if we do not believe we know any of 
the three terms (X-ray luminosity, radio luminosity and BH 
mass) as being the primary physical driver, all should be 
treated equal in a correlation study (but this is beyond the 
purposes of this paper). 

In order to improve these analysis it would be very use- 
ful to obtain deeper radio observations of complete samples 
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of AGNs combined with detailed optical-NIR-MIR SED ob- 
servations. These data, when available, will eventually allow 
to measure the dependence of the radio luminosity distribu- 
tion (i.e. the feedback) from the star and BH masses, their 
derivatives (star formation and accretion rates) and the red- 
shift. A result that could be achieved by complementing mul- 
tiwavelength surveys with observations carried out with new 
radio facilities such as the Expanded Very Large Array and 
the Square Kilometer Array. 
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